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ABSTRACT 

The inner mass-profile of the relaxed cluster Abell 1703 is analysed by two very 
different strong-lensing techniques applied to deep ACS and WFC3 imaging. Our 
parametric method has the accuracy required to reproduce the many sets of multiple 
images, based on the assumption that mass approximately traces light. We test this 
assumption with a fully non-parametric, adaptive grid method, with no knowledge of 
the galaxy distribution. Differences between the methods are seen on fine scales due 
to member galaxies which must be included in models designed to search for lensed 
images, but on the larger scale the general distribution of dark matter is in good agree- 
ment, with very similar radial mass profiles. We add undiluted weak- lensing measure- 
ments from deep multi-colour Subaru imaging to obtain a fully model-independent 
mass profile out to the virial radius and beyond. Consistency is found in the region 
of overlap between the weak and strong lensing, and the full mass profile is well- 
described by an NFW model of a concentration parameter, c v i r c± 7.15 ± 0.5 (and 
M vir ~ 1.22 ± 0.15 x 10 15 M Q /h). Abell 1703 lies above the standard c-M relation 
predicted for the standard ACDM model, similar to other massive relaxed clusters 
with accurately determined lensing-based profiles. 

Key words: dark matter, galaxies: clusters: individuals: Abell 1703, galaxies: clus- 
ters: general, gravitational lensing , galaxies: elliptical and lenticular, cD, galaxies: 
formation 



1 INTRODUCTION 

Simulated CDM dominated halos consistently predict mass 
profiles that steepen with radius, providing a distinctive, 
fundamental prediction for this form of dark matter (DM; 
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Navarro, Frenk & White 1996; NFW). Furthermore, the de- 
gree of mass concentration should decline with increasing 
cluster mass because in the hierarchical model massive clus- 
ters collapse later, when the cosmological background den- 
sity is lower. These predictions are now being subject to 
stringent lensing based analyses, using multiply- lensed im- 
ages and with weak lensing (WL) information. To date only 
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a few clusters have been reliably analysed by combining both 
weak and strong lensing for a full determination of the mass 
profile and a definitive comparison with predictions, (e.g., 
Gavazzi et al. 2003, Broadhurst et al. 2005b, 2008, Umetsu 
et al. 2010, Newman et al. 2009, Okabe et al. 2010 and ref- 
erences therein). The upcoming multi-cycle HST program 
of cluster imaging (the CLASH program 1 ) will provide a 
much more definitive derivation of mass profiles for a statis- 
tical sample of relaxed, X-ray selected clusters, combining 
high resolution space imaging with deep, wide- field ground 
based data. 

Strong gravitational lensing (SL) is of great signifi- 
cance as a cosmological probe, providing model-free masses 
of galaxies and clusters interior to the Einstein radius and 
useful constraints on their inner mass profiles. The mass 
density in the central regions of distant clusters typically 
exceeds the critical value required for lensing, generating 
multiple-images of background objects (e.g., Horesh et al. 
2010, Kausch et al. 2010). Recent analyses have shown that 
many sets of multiply-lensed images can be uncovered with 
high-quality space imaging and thanks to improved mod- 
elling techniques. Reliable mass maps are claimed for sev- 
eral well studied clusters with deep space imaging (e.g., 
Abell 370, Kneib et al. 1993, Richard et al. 2010; Abell 
611 Newman et al. 2009; Abell 901, Deb et al. 2010; Abell 
1689, Broadhurst et al. 2005a, Coe et al. 2010; C10024+1654, 
Liesenborgs et al. 2008, Zitrin et al. 2009b; MS 2137.3-2353, 
Gavazzi et al. 2003, Merten et al. 2009; RXJ1347, Bradac et 
al. 2008, Halkola et al. 2008; SDSS J1004+4112, Sharon et 
al. 2005; "The bullet cluster", Bradac et al. 2006). 

It is important to realise that most published mass maps 
usually adopt an initial model gradient for the cluster mass 
profile, rather than deducing and constraining it directly 
from the data. Various SL modelling methods have devel- 
oped over the past two decades in response to the huge im- 
provements in astronomical imaging. Most methods can be 
classified as "parametric" if based on model prescriptions, 
or "non-parametric" if "grid-based" or interpolative, capable 
of arbitrary forms (e.g., Saha & Williams 1997, Abdelsalam, 
Saha, & Williams 1998, Diego et al. 2005, Liesenborgs et al. 
2006, Valls-Gabaud et al. 2006; see also §4.4 in Coe et al. 
2008). The non-parametric grid methods do not have the 
resolution to accurately locate and reproduce multiple im- 
ages, and usually rely on images identified by other means 
- often just eyeball candidates or those identified from the 
subset of parametric models with predictive power to locate 
images, for which the number of free parameters does not 
exceed the number of independent multiple-images used as 
constraints. 

The method developed by Broadhurst et al. (2005a), 
and simplified further by Zitrin et al. (2009b) has securely 
identified tens of multiple images in high quality ACS im- 
ages, behind Abell 1689 and C10024 and also a sample of 
12 MACS clusters at z > 0.5 sample (Zitrin et al. 2010), 
with only six free parameters so that in practice the number 
of multiple images uncovered readily exceeds the number of 
free parameters, as minimally required in order to get a re- 
liable fit. This approach to SL is based on the assumption 
that mass approximately traces light, and will be employed 
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here as our parametric analysis of Abell 1703. We also ap- 
ply the non-parametric technique of Liesenborgs et al. (2006, 
2007, 2009) which employs an adaptive grid inversion tech- 
nique and has been well tested on available multiple-imaging 
data such at the many sets of multiple-images uncovered by 
Zitrin et al. (2009b; with photometric redshifts calculated 
therein), in C10024+1654 (see also Liesenborgs et al. 2008). 

Here we compare in more detail these two very differ- 
ent methods applied to Abell 1703, allowing in principle a 
test of the assumptions behind the parametric technique. 
Abell 1703 has been subject to various complementary stud- 
ies ranging from early X-ray and optical work to more recent 
high quality lensing analyses from ground and space (Kowal- 
ski et al. 1984, Bohringer et al. 2000, Leir & Van Den Bergh 
1977, Bade et al. 1998, Cooray et al. 1998, Koester et al. 
2007, Scott et al. 2007, Riley et al. 2008,2009, Bruursema 
et al. 2008, Oguri et al. 2009), with detailed radio sources 
(Rizza et al. 2003, Coble et al. 2007) and highly magnified 
high- 2 galaxies (Zheng et al. 2009). SL analyses of Abell 1703 
have been carried out by Hennawi et al. (2008), Limousin et 
al. (2008), de Xivry & Marshall (2009), Oguri et al. (2009), 
Richard et al. (2009), and Saha & Read (2009). 

The cluster Abell 1703 (Abell 1958; see also Abell, Cor- 
win, & Olowin 1989) is known to have many sets of multiple 
images with impressive spectroscopic redshift information 
(Limousin et al. 2008, Richard et al. 2009; see also Estrada 
et al. 2007, Hennawi et al. 2008). Here we take advantage 
of two independent SL modelling techniques which are in- 
teresting to compare given their very different approaches. 
We then add the accurate WL data from Broadhurst et al. 
(2008; see also Medezinski et al. 2010) to complete the mass 
profile for comparison with theoretical predictions out to the 
virial radius and beyond and to examine the consistency of 
the WL and SL derived profiles in the region of overlap. 

The paper is organised as follows: In §2 we describe 
the observations. In §3 we detail the SL modelling meth- 
ods and their implementation. In §4 we report and discuss 
the results, which are then summarised in §5. Throughout 
this paper we adopt a concordance ACDM cosmology with 
(Q m0 = 0.3, Q A0 = 0.7, h = 0.7). We adopt a redshift of 
z = 0.28 for the cluster, equal to that of the prominent 
central BCG galaxy (Allen et al. 1992). With these param- 
eters one arcsecond corresponds to a physical scale of 4.25 
kpc for this cluster. The reference centre of our analysis is 
fixed at the centre of the BCG: RA = 13:15:05.24, Dec = 
+51:49:02.6 (J2000.0). 



2 OBSERVATIONS 

Abell 1703 was observed in November 2004, with the Wide 
Field Channel (WFC) of the ACS installed on HST, in 
the framework of the ACS Guaranteed Time Observations 
(GTO; Ford et al. 2003) which includes deep observations 
of several massive, intermediate-redshift galaxy clusters. In- 
tegration times of 7050, 5564, 5564, 8494+1340, 5564 x 2 
and 8900 x 2 seconds, were obtained through the F435W, 
F475W, F555W, F625W, F775W, and F850LP filters, re- 
spectively, and are available in the Hubble Legacy Archive. 
Some important aims of the GTO program are determina- 
tion of the mass distribution of clusters for testing the stan- 
dard cosmological model and to study distant, background 
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Figure 1. The starting point of the parametric model, where we 
define the surface mass distribution based on the cluster member 
galaxies (see ^3} . Axes are in ACS pixels (0.05" /pixel). North is 
up, East is left. 



lensed galaxies for which some of the very highest redshift 
galaxies are known because of high magnification by massive 
clusters (Franx et al. 1997, Frye & Broadhurst 1998, Frye, 
Broadhurst & Bem'tez 2002, Kneib et al. 2004, Stark et al. 
2007, Bouwens et al. 2009, Bradley et al. 2008, in Abell 1703: 
Zheng et al. 2009). 

As part of the ACS GTO cluster program, in April 2010 
we also observed Abell 1703 with the near-infrared channel 
of HST new Wide Field Camera 3 (WFC3/IR). The obser- 
vations consisted of 1 orbit (2812s) each in the F125W and 
F160W bands. 

Various redshifts have been quoted for Abell 1703 (e.g., 
Struble & Rood 1987,1999, Bohringer et al. 2000), corre- 
sponding to several different cluster galaxies, out of which 
we adopt that of the prominent BCG galaxy at z = 0.28 
(Allen et al. 1992). This redshift was also used in recent SL 
work on this cluster by Limousin et al. (2008) and Richard et 
al. (2009) who identified many sets of multiple-images which 
we incorporate in this work, as will be detailed below. 



3 STRONG LENSING MODELLING AND 
ANALYSES 

3.1 Parametric Method 

We first apply our well tested approach to lens modelling, 
which has previously uncovered large numbers of multiply- 
lensed galaxies in ACS images of Abell 1689, C10024, and 
12 high-z MACS clusters (respectively, Broadhurst et al. 
2005a, Zitrin et al. 2009b, Zitrin & Broadhurst 2009, Zitrin 
et al. 2009a, 2010). The full details of this approach can be 
found in these earlier papers. Briefly, the basic assumption 
adopted is that mass approximately traces light, so that the 
photometry of the red cluster member galaxies is used as 
the starting point for our model. Cluster member galaxies 
are identified as lying close to the cluster sequence by the 
photometry provided in the Hubble Legacy Archive. 

We approximate the large scale distribution of cluster 
mass by assigning a power-law mass profile to each galaxy 



Figure 2. The resulting smooth mass component of the paramet- 
ric model (see ^3}. Axes are in ACS pixels (0.05 /; /pixel). North 
is up, East is left. 

(see Figure [T]), the sum of which is then smoothed (see Fig- 
ure [5]). The degree of smoothing (S) and the index of the 
power-law (q) are the most important free parameters de- 
termining the mass profile. A worthwhile improvement in 
fitting the location of the lensed images is generally found 
by expanding to first order the gravitational potential of 
this smooth component, equivalent to a coherent shear de- 
scribing the overall matter ellipticity, where the direction of 
the shear and its amplitude are free parameters, allowing 
for some flexibility in the relation between the distribution 
of DM and the distribution of galaxies, which cannot be ex- 
pected to trace each other in detail. The total deflection field 
ar(0), consists of the galaxy component, a ga i(0), scaled by 
a factor K ga i, the cluster DM component c*dm(#), scaled by 
(1-Kg a i), and the external shear component a ex (0): 

a T (0) = K ga ia ga i(6) + (1 - K ga i)a DM ($) + a ex (0), (1) 

where the deflection field at position m due to the external 
shear, a ex (6 m ) = (a eXjX , a ex , y ), is given by: 

a ex ,x(0m) = |tI cos(20 7 )Ax m + M sin(20 7 )Ay m , (2) 

aex, y (0 m ) = |7| sin(20 7 )Ax m - | -y | cos(2</> 7 ) Ay m , (3) 

where (Ax m , Ay m ) is the displacement vector of the posi- 
tion m with respect to a fiducial reference position, which 
we take as the lower- left pixel position (1,1), and 7 is the 
position angle of the spin- 2 external gravitational shear mea- 
sured ant i- clockwise from the x-axis. The normalisation of 
the model and the relative scaling of the smooth DM compo- 
nent versus the galaxy contribution brings the total number 
of free parameters in the model to 6. This approach to SL 
is sufficient to accurately predict the locations and internal 
structure of multiple images, since in practice the number 
of multiple images uncovered readily exceeds the number of 
free parameters thus fully constraining them. 

In addition, two of the six free parameters can be pri- 
marily set to reasonable values so only 4 of these parame- 
ters have to be constrained initially, which sets a very re- 
liable starting-point using obvious systems. The mass dis- 
tribution is therefore primarily well constrained, uncovering 
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Figure 3. Galaxy cluster Abell 1703 [z = 0.28) imaged with Hubble/ACS. The overlaid critical curve (blue) corresponds to the distance 
of system 1, at z s = 0.889. The outer white critical curve corresponds to the giant arc (systems 10/11) at higher redshift, z s = 2.627, 
enclosing a critical area of an effective Einstein radius of ~ 130 kpc at the redshift of this cluster. North is up, East is left. 



many multiple-images which can be then iteratively incor- 
porated into the model, by using their redshift estimation 
and location in the image-plane. 

Firstly we use this preliminary model to lens the more 
obvious lensed galaxies back to the source plane by subtract- 
ing the derived deflection field, and then relens the source 
plane to predict the detailed appearance and location of ad- 
ditional counter images, which may then be identified in 
the data by morphology, internal structure and colour. We 
stress that multiple images found this way must be accu- 
rately reproduced by our model and are not simply eyeball 
"candidates" requiring redshift verification. In Abell 1703 
many multiple-images (16 systems, most of them with spec- 
troscopic redshifts; Limousin et al. 2008, Richard et al. 2009) 
are already known and therefore simply used to constrain 
the fit, which is assessed by the RMS uncertainty in the 
image plane: 

RMSi mages = y^((gj - Xi) 2 + (y % - yif) / N images , (4) 

where x i and y i are the locations given by the model, and 
Xi and yi are the real images location, and the sum is over 



all Nimages images. The best-fit solution is unique in this 
context, and the model uncertainty is determined by the 
location of predicted images in the image plane. Impor- 
tantly, this image-plane minimisation does not suffer from 
the well known bias involved with source plane minimisa- 
tion, where solutions are biased by minimal scatter towards 
shallow mass profiles with correspondingly higher magnifi- 
cation. 

The model is successively refined as additional sets of 
multiple images are incorporated to improve the fit, using 
also their redshift measurements or estimates, for better con- 
straining the mass slope through the cosmological relation 
of the Dis/Ds growth. 

3.2 Non-parametric inversion method 

3.2.1 Genetic algorithm based inversion 

The non-parametric inversion method that we apply here is 
based on the work of Liesenborgs et al. (2006). It requires the 
user to specify a square-shaped region in which the inversion 
routine should try to reconstruct the mass distribution. Ad- 
ditionally, it is necessary to define which images correspond 
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Figure 4. Distribution of multiply- lensed images used to constrain the models marked on a colourful Hubble/ ACS image of Abell 1703, 
with some central galaxies including the cD galaxy subtracted, and its original location marked in red. Subtraction has been carried 
out by modelling the main cD galaxy and three other galaxies using Chebyshev- Fourier basis functions ( "CheF-lets" ; Jimenez-Teja Sz 
Benftez, in preparation). These multiply-lensed images and more details can be found in Limousin et al. (2008) and Richard et al. (2009). 
North is up, East is left. 



to the same source and at what redshifts the sources are 
located. In a first step, the square region is subdivided in a 
uniform way into a number of smaller square grid cells, and 
to each cell a projected Plummer sphere (Plummer, 1911) 
is assigned. The width of each basis function is set propor- 
tional to the grid cell size. As an additional basis function, 
a sheet of mass can be included; this can be useful as in 
the center of clusters a non-negligible density offset may be 
present which can prove difficult to model using Plummer 
basis functions. A genetic algorithm is then used to search 
for appropriate weights of these basis functions, yielding a 
first approximation of the projected density of the lens. 

Using this first approximate solution, a new grid is then 
constructed in which regions containing more mass are sub- 
divided further. It should be noted that the mass sheet ba- 
sis function is not taken into consideration in this step as 
it is structureless. Using this new grid, basis functions are 
assigned and the genetic algorithm again looks for appro- 
priate weights. This refinement procedure can be repeated 
until the added resolution no longer results in an improved 
fit to the data. 

The actual search for appropriate weights of the basis 
functions and thus for the mass distribution, employs a ge- 
netic algorithm. This is a heuristic optimisation strategy, 
inspired by the theory of evolution by Darwin. In essence, 



one tries to breed solutions to a problem, by evolving an 
initial population of trial solutions towards solutions which 
are better adapted to the problem under study. To create 
the next generation from the current one, trial solutions are 
combined and mutated, while applying selection pressure, 
i.e., making sure that solutions which are deemed better 
create more offspring. In this approach it is even possible 
to simultaneously optimise against several so-called fitness 
measures; one then speaks of a multi-objective genetic algo- 
rithm (see, e.g., Deb 2001). 



As the number of basis functions used can become quite 
large (e.g. ~ 1000) and the genetic algorithm starts from 
random initial solutions, several runs of this procedure will 
produce results that differ somewhat. Therefore a set of so- 
lutions is usually generated, so that one can inspect the 
common features of these mass maps, and the standard de- 
viation can be used as a measure of the reliability at any 
location. The algorithm details and original fitness criteria 
are described in Liesenborgs et al. (2006, 2007, 2009). Be- 
low, the fitness criteria used in this work shall be described. 
It is based on these criteria that selection pressure will be 
applied. 
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counted and the total count for all sources is used as the null 
space fitness measure. As for each source this gives an ap- 
proximation of the amount of images, a lower value indicates 
a better fitness with respect to this criterion. 



3.2.3 Input 

The input of the inversion routine consists of multiply im- 
aged systems together with their redshifts. The systems 
listed in Richard et al. (2009) were used for this purpose, 
and where available, spectroscopic redshift information was 
used. For systems 2, 8, 9 and 12, the redshift predicted by 
the model in this work was used, as these redshifts were in 
good agreement with the photometric ones. Note, systems 
13 and 14 were not used in the inversion as the redshift 
estimates seemed more uncertain (however we do find that 
including them results only in minor changes to the mass 
model). Based on these image systems, the inversion rou- 
tine was instructed to look for mass in a 2x2 arcmin 2 region, 
roughly centered on the cD galaxy. To limit the amount of 
predicted images that were not part of the input, and there- 
fore to limit the amount of unnecessary substructure, the 
null-space region was 3x3 arcmin 2 in size, and a 48x48 grid 
was used for calculating the corresponding fitness measure. 



Figure 5. An optical/ ACS + IR/WFC3 colourful image of the 
central region of Abell 1703. The much wider colour range enabled 
by incorporation of the IR data manifests the distinct colours of 
the many multiple-systems shown in Figures [3] and ^] and the 
identification of an additional system (17) marked in circles, fol- 
lowing the same symmetry of systems 15 and 16 (Figure^}. 

3.2.2 Fitness criteria 

In strong lens inversion one tries to deduce the projected 
mass distributions based on data of multiply imaged sources. 
Since each set of images originates from a single source, 
projecting the images back onto the corresponding source 
plane should produce a consistent source. Previously, only 
extended images could be used, in which case the back- 
projected images should overlap in the source plane. To 
calculate the amount of overlap, the estimated size of the 
source was used as a length scale. The method was adapted 
to work with point images as well. In this case the envelope 
of all estimated source positions is used as a length scale 
when calculating the distances between the back-projected 
images of each source. Using this length scale instead of an 
absolute scale, prevents scaling the source plane to obtain 
a better fitness value. This is especially important when a 
mass sheet is included as a basis function since it has pre- 
cisely this effect. 

Apart from the locations where images can be seen, ad- 
ditional constraints come from the area in which no images 
are observed, i.e., the null space. To avoid predicting ex- 
tra images, which corresponds to avoiding unnecessary sub- 
structure, the user can define a region which will be used 
to check for additional images. This region is divided into 
a large number of triangles and for each source, the trian- 
gles are projected onto the source plane. Each triangle that 
overlaps with the envelope of the back-projected images is 



4 RESULTS AND DISCUSSION 

4.1 Strong-Lensing Regime and cD Galaxy 

In the SL regime we have modelled Abell 1703 (see Fig- 
ure [3]) using the many sets of multiply- lensed images previ- 
ously identified, most of which have spectroscopic redshifts 
reported in Limousin (2008) and Richard et al. (2009). Here 
we incorporate also new WFC3/IR imaging, which reveals 
the distinct colours of each system, and enables the identifi- 
cation of an additional system (17) following the same sym- 
metry as systems 15/16 (see Figures [4] and [5]) with a similar 
model redshift of z s ~ 2.8. We have used the parametric 
method of Zitrin et al. (2009b) to further verify the reliabil- 
ity of the many multiply-lensed images across the field and 
to securely input them into the non-parametric method of 
Liesenborgs et al. (2006) to model the central mass distri- 
bution. 

Our parametric model (see Figure [6]) accurately repro- 
duces all multiply-lensed images, indicating that our pre- 
liminary assumption that mass traces light is reasonable. In 
order to further test this assumption we then applied the 
non-parametric technique of Liesenborgs et al. (2006, 2007, 
2009) for which no prior information for the distribution of 
cluster galaxies or mass is input. Nevertheless, the results 
of this method seem to trace the distribution of light as 
can be seen in Figure [7] generating a 2D mass distribution 
which is remarkably similar to the result of the parametric 
mass model on the large scale, shown in Figure [6] A distinct 
substructure is seen in both maps and corresponds to local 
galaxy over-densities. In addition, the two methods produce 
very similar mass profiles over a range of scales covering the 
full distribution of multiple images, with a mean inner slope 
of dlog£/dlog6> ~ -0.5. 

We have examined the difference between these two 
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Figure 6. 2D surface mass distribution («), in units of the critical 
density (for z s = 2.627), of Abell 1703. Contours are shown in 
linear units, derived from the parametric mass model constrained 
using the many sets of multiply-lensed images seen in Figure ^] 
Axes are in ACS pixels (0.05" /pixel), North is up, East is left. 

mass maps by subtracting the non-parametric mass distri- 
bution from the parametric mass distribution. The result 
is shown in Figure [8] As can be seen in this figure, the 
main positive differences, marked in red, occur mainly where 
galaxies are located in the data, since these must contain 
mass and are included only in the parametric model. The 
non-parametric model does place mass at these locations, 
but it is usually smoothly distributed as this approach to 
modelling does not make prior assumptions about the mass 
distribution and thus does not achieve a spatial resolution 
sufficient for resolving individual cluster galaxies. The main 
negative differences, marked in blue, are seen where the non- 
parametric model has more mass than implied by the galaxy 
distribution, but overall these are small and likely inevitable 
given the inherent noise set by the finite amount of input 
data. The mean difference across this field is \Ak\ — 0.19, 
contributed mainly by the inclusion of cluster members or 
discrepancies outside the critical curves, where one has rel- 
atively poor constraints from the observed multiple images. 

The parametric method of Zitrin et al. (2009b) has been 
shown to have the predictive power to find many multiple 
images in the field. This parametric method has inherently 
more structure on small scales by virtue of the inclusion 
of cluster members which can significantly deflect images 
locally and must be included in order to find lensed im- 
ages. Due to the low number of parameters this model is 
initially well-constrained using only a few sets of usually ob- 
vious multiple-systems, thus correlated to the initial mass 
distribution so that the image-plane reproduction accuracy 
can be only somewhat improved as newly-found multiple- 
systems are incorporated, but the overall gradient of the 
cluster lensing profile is significantly refined through the cos- 
mological distance-redshift relation. It is important to have 
a wide range of background source redshifts for a reliable 
profile determination, otherwise the SL models are degener- 
ate with respect to the profile, although the relative distri- 




Figure 7. 2D surface mass distribution (ac) contours overlayed on 
the cluster image, derived from the non-parametric mass model 
constrained using the multiply-lensed images seen in Figure^] As 
can be seen, though no galaxies were included in this modelling 
method, the mass contours steepen up where significant galaxies 
are present. Clearly this mass distribution is in good agreement 
with the parametric mass distribution seen in Figure [6] 

bution of matter and substructure can still be reasonable, 
finding many sets of multiple-images. In the non-parametric 
approach on the other hand, the fit is much more flexible and 
is continuously improved as more of these images are incor- 
porated and the overall solution is clarified, allowing the ex- 
clusion of a wide range of non-unique solutions (Liesenborgs 
et al. 2006, 2008). When sufficient images are incorporated, 
the overall mass distributions and profiles of these two meth- 
ods become very similar, as we have found here and is shown 
in the comparison of Figures [6] |7l and [8] 

The critical curves for different sources are plotted on 
the cluster image in Figure [3] The critical curves for a source 
redshift of z s = 2.627 enclose an area with an effective Ein- 
stein radius of 30.5 =b 3" (—130 kpc at the redshift of the 
cluster), and a mass of 1.25 ± 0.1 x 10 14 M o . Our paramet- 
ric model reproduces all multiply-lensed systems within a 
Gi — 1.5" of their real location. In particular, our best- 
fit parametric model accurately reproduces the complicated 
ring (z s — 0.889; system 1) as can be seen in Figure [9] We 
notice that only models with central mass-profiles steeper 
than a certain threshold accurately reproduce all parts of 
the ring, importantly enabling us to constrain the mass and 
the profile of the cD galaxy in the range c± 1 — 5", where the 
closest image of the ring system forms. This is shown also 
in Figure [To] 

We find that the cD galaxy encloses a projected mass 
of 5.2 ± 0.4 x 10 n M o within a radius of ~ 5" (-22 kpc) 
after subtracting the interpolated smooth DM component 
(— 6.3 x 1O 12 M inside this aperture), and has a B-band 
luminousity of 8.8 ± 0.1 x 1O 1O L (fluxes were converted to 
luminousities using the LRG template described in Bemtez 
et al. 2009). This corresponds to an averaged M/Lb of 
~ 6 (M/L)q in this region. This ratio can be fully accounted 
for by the stars contained in this galaxy, for which we obtain 
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Figure 8. A 2D contour map of the difference between the para- 
metric (Figure |SJ and the non-parametric (Figure [7| mass dis- 
tributions. Contours are in A«, plotted in equal linear spaces of 
0.1. The main positive differences (red) are seen where galaxies 
are located in the data (since these are included only in the para- 
metric model), and the main negative differences (blue) are seen 
where the non-parametric model has more mass than implied by 
the galaxy distribution. Still these differences are overall small, 
with a mean difference of | Aft | =0.19 across this field, contributed 
mainly by the inclusion of cluster members or discrepancies out- 
side the critical curves, where one has relatively poor constraints 
from the observed multiple images. 



as well M/Lb — 6 (M/L) , for a single-burst stellar popu- 
lation formed at z = 3 and viewed at a redshift of z = 0.28, 
equivalent to an age of ~8.1 Gyrs, and with half solar metal- 
licity (by evolutionary models of Bruzual & Chariot 2003). 
This result is in agreement with the result of Limousin et 
al. (2008), and is similar to other lensing based cD masses 
in well-studied clusters, for which low M/L ratios are also 
found and fully accounted for by the measured stellar light 
(e.g., Gavazzi et al. 2003 for MS2137-2353, and Zitrin & 
Broadhurst 2009 for MACS J1149. 5+2223). 

We mentioned in the preceding sections that the pro- 
file can only be accurately constrained by incorporating the 
cosmological redshift-distance relation, i.e., the lensing dis- 
tance of each system based on the measured spectroscopic 
redshifts. In so doing we make use in particular of the 
z — 0.889 system (the ring; system number 1), whose red- 
shift is very distinct from the rest of the multiple-image 
systems, thus strongly constraining the profile. We exam- 
ine how well the cosmological relation is reproduced by the 
parametric model, accounting also for all other systems with 
reliable spec-z measurements, as shown in Figure [TT] Clearly 
the redshifts of these systems verify very well that the pre- 
dicted deflection of the best fitting model at the redshift of 
each of these systems, lies precisely along the expected cos- 
mological relation, with a mean deviation of only Af < 0.01 
(see Figure 11), and % 2 — 0-1 f° r the best model, consid- 





Figure 9. Reproduction of system 1 by our model, by delensing 
image 1.4 into the source plane, and then relensing the source 
plane pixels onto the image plane to accurately form the ring. 
By tuning the inner slope of the mass distribution, the observed 
structure is reproduced very closely. A small fifth image, 1.5, is 
formed outside the tangential critical curve (for z s = 0.889; see 
Figures [3]and[4]). 
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erably strengthening the plausibility of our parametric ap- 
proach to modelling in general. In Figs |12| [l3l and ^] we 
give further examples demonstrating how different systems 



Figure 10. The four images of the ring (system 1, z s = 0.889; 
see also Figure [9} next to the cluster core enable a unique de- 
termination of the cD galaxy mass-profile in the range ~ 1 — 5 /x 
(i.e., ~ 4 — 22 kpc; where the closest image of the ring appears). 
The solid thick black line represents the best-fit parametric model 
which accurately reproduces the ring. The dashed red line repre- 
sents the minimum "shallowness" threshold, meaning that only 
models higher than this threshold will form the four images of the 
ring, and the dash-dotted blue line shows the limit above which 
these images become overly distorted. 



are accurately reproduced by our parametric model, in addi- 
tion to the remarkable reproduction of the ring system seen 
in Fig. [9] 

4.2 Combined Weak and Strong Lensing 

We now compare our SL analysis results with undiluted 
WL data from deep Subaru g'r'i' images (Broadhurst et 
al. 2008). The Subaru WL data, covering a wide field of 
« 34' x 27', allow us to probe the cluster mass distribution 
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Figure 11. Growth of the scaling factor f(di s /d s ) as a function 
of redshift, normalised so f=l at z = 2.627. Plotted lines are 
the expected ratio from the chosen specified cosmological model. 
The circles correspond to the multiple-image systems reproduced 
by the parametric mass model, versus their real spectroscopic 
redshift. The data follow very well the relation predicted by the 
standard cosmological model (mean deviation of only Af < 0.01, 
and x 2 = 0.1 for this fit). 




Figure 13. Reproduction of system 6 by our model, by delensing 
image 6.2 into the source plane, and relensing the source-plane 
pixels onto the image plane to accurately form the resolved inter- 
nal details of the other images of this system. 





Figure 14. Reproduction of the giant arc, systems 10 Sz 11, by 
our model. We delens image 10.2 into the source plane, and relens 
the source-plane pixels onto the image plane, which accurately 
reproduces also the other half of the main arc, image 10.1, and 
two other smaller images (10.3 and 10.4; see Figure[3}. 



Figure 12. Reproduction of systems 4 Sz 5 by our model, by 
delensing image 4.1/5.1 into the source plane, and relensing the 
source-plane pixels onto the image plane to accurately form the 
other images of this system. 



over a wide radial range, 6 ~ [0.7', 18'], in the subcritical 
regime (0 > 30"). 

For a direct comparison with the WL data, we followed 
the method outlined in Umetsu et al. (2010) to translate 
our SL mass profiles into corresponding tangential distortion 
profiles g+(0) = 7+(0)/(l — k(0)) for a fiducial source red- 
shift z s = 1, roughly matching the mean depth of blue+red 
background galaxies used for the WL analysis. In Figure 
|15| we compare our parametric and non-parametric SL in- 
ner profiles with the Subaru distortion profile. Our SL and 
WL results are in good agreement where the data overlap, 
= [40", 90"], for both SL methods. Furthermore, a simple 
inward extrapolation of the best-fitting NFW profile (see 
Table [T]) for the outer Subaru observations with input of 
the Einstein radius, fits well with the inner SL information, 
in particular, slightly better for the non-parametric profile. 
The parametric profile has a minor "bump" around ~ 30" 



due to other bright galaxies in the field that are not in- 
cluded in the non-parametric model. This translates into a 
dip in the g+(0) profile of the parametric technique, slightly 
deviating from the smooth NFW curve. Furthermore, the 
WL+Einstein radius NFW fitting assumes a circularly sym- 
metric lens, possibly biasing this NFW profile fit done this 
way. 

We then reconstruct the outer mass profile from the 
Subaru WL data using the shear-based one-dimensional in- 
version method outlined in Umetsu et al. (2009, 2010). Fig- 
ure [16] compares our SL and WL results in terms of the lens 
convergence profile, ft(0), where the combined SL and WL 
results produce a coherent mass profile with a continuously 
steepening radial trend from the central region to the out- 
skirts of the cluster. 

Overall, we combine our SL and WL results to examine 
the form of the underlying cluster mass profile and to char- 
acterise cluster mass and structure properties. To do this, we 
fit our SL and WL constraints with an NFW profile in four 
independent manners: First, we fit the Subaru distortion 
profile g+(0) alone with no SL information involved. Second, 
we utilise the inner Einstein radius constraint, $e = 30. 5 ±3" 
(10% uncertainty) at z s = 2.627, in conjunction with the 



10 Zitrin et al 



Subaru distortion profile (for details, see Umetsu & Broad- 
hurst 2008 and Umetsu et al. 2010). Most credibly, we then 
fit our joint SL+WL convergence profiles k(0) of the para- 
metric and non-parametric SL methods following the pre- 
scription given by Umetsu et al. (2010). For the joint NFW 
fitting we constrain the SL data to the range [5", 25"] where 
the two independent SL profiles are extremely similar, and 
since in general the shapes of the SL mass profiles within the 
Einstein radius look smoother, matching better the expec- 
tation from the outer WL profile. All these methods yield 
similar and consistent results, which are also summarised in 
Tabled] 

The NFW universal density profile has a two-parameter 
functional form as (Navarro, Frenk, & White 1997) 



r [h *kpc] 

200 500 1 ooo 



PNFw(r) 



(r/r s )(l+r/r s ) 2 



(5) 



where p s is a characteristic inner density, and r s is a 
characteristic inner radius. The logarithmic gradient n = 
d In p(r) /d In r of the NFW density profile flattens continu- 
ously towards the center of mass, with a flatter central slope 
n = — 1 and a steeper outer slope (n — ► —3 when r — ► oo) 
than a purely isothermal structure (n = —2). A useful in- 
dex, the concentration, compares the virial radius, r v i r , to 
r s of the NFW profile, c v i r = r v i T /r a . 

We specify the NFW model with the halo virial mass 
M v ir and the concentration c v i r instead of p s and r s . Here 
the errors for these best-fit NFW parameters include the 
uncertainty in the source redshift calibration for WL, z s = 
1.0 ±0.2. The typical halo concentration obtained by the four 
methods described above IS Cvir ~ 7.5 =b 0.5, and the typical 
virial mass is M vir ~ 1.2 ± 0.15 x 1O 15 M h' 1 (r v i r ~ 1.84 
Mpc/i" 1 ). The joint SL+WL fits yield c vir ~ 7.15 ± 0.5 
and M vir ~ 1.22 ± 0.15 x 10 15 M o /i _1 , agreeing with the 
recent result of Oguri et al. (2009; c v ir — 6.51J" 7, M V i r ~ 
1.05lo.25 x lO^Mo/i -1 ). This puts Abell 1703 above the 
standard c-M relation and manifests again the tension with 
the standard ACDM model. This can be seen in Figure |T7] 
where we plot confidence levels of the concentration param- 
eter derived for Abell 1703, along with c-M relations includ- 
ing la uncertainties, deduced from simulations by Duffy et 
al. (2008) with WMAP5 parameters and scaled to z c = 0.28. 
Since the concentration parameter depends on cluster forma- 
tion time (as discussed recently by Sadeh & Rephaeli 2008) 
and though such conclusions have previously been reached 
for several other massive clusters (e.g., Broadhurst et al. 
2008, Umetsu et al. 2010), a similar comparison should nat- 
urally be made for many other clusters in order to clearly 
establish that a statistically significant trend is discerned. 



5 SUMMARY 

In this work we have produced an accurate full-range profile 
of the cluster Abell 1703, based on combined WL and SL 
mass models applied to exceptionally high quality space and 
ground based imaging data, including recent observations of 
the WFC3/IR on HST. This profile adds to the relatively 
few clusters for which precise and reliable mass profile have 
been constructed, including MS2137 (Gavazzi et al. 2003, 
Merten et al. 2009), Abell 1689 (Broadhurst et al. 2005a,b, 
Lemze et al. 2008, Umetsu & Broadhurst 2008, Limousin et 
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Figure 15. Comparison of our parametric and non-parametric SL 
profiles with the Subaru distortion profile. Our SL and WL results 
are in good agreement where the data overlap, = [40", 90"], for 
both SL methods. Furthermore, a simple inward extrapolation of 
the best-fitting NFW profile (black solid curve; see also Table [TJ 
for the outer Subaru observations with input of the Einstein ra- 
dius, fits well with the inner SL information, in particular, slightly 
better for the non-parametric profile. The parametric profile has 
a minor "bump" around ~ 30" due to other bright galaxies in the 
field which are not included in the non-parametric model, trans- 
lating into a dip in the g+ (6) profile of the parametric technique, 
slightly deviating from the smooth NFW curve. 
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Figure 16. Comparison of our SL and WL results in terms of the 
lens convergence profile, k{9). The combined SL and WL results 
produce a coherent mass profile with a continuously steepening 
radial trend from the central region to the outskirts of the cluster. 
Both SL methods are in very good agreement with the WL data 
and the NFW fit (see also Table [TJ. 



al. 2007), Abell 2218 (Kneib et al. 1996, Abdelsalam, Saha 
& Williams 1998) and C10024 (Zitrin et al. 2009b, Umetsu 
et al. 2010). 

We have applied two independent SL techniques in or- 
der to derive the inner mass profile and to examine for con- 
sistency the basic assumption of the parametric-model, that 
mass traces light. Both techniques derive remarkably sim- 
ilar mass distributions with the same major substructure 
present, and with similar overall radial mass profiles (with 
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Table 1. Comparison of the different NFW parameters obtained independently by different lensing techniques. Column 1: The method; 
Column 2: Resulting concentration parameter, c v i r ; Column 3: Resulting virial mass, M v i r , in 1O 15 M0//i; Column 4- Reduced x 2 of 
the fit; Column 5: Q- value goodness of fit. 



Method 


Cyir 


M v i r 


Red. x 2 


Q-value 


WL ± Parametric SL 
WL ± Non-parametric SL 
WL alone 

WL + e e 


7.07 ±0.47 
7.23 ±0.45 

7 71 +2.22 
'•' 1 -1.66 
7g2 +1.34 


-1 ,o+0.12 

1 - 18 -o.n 

1 26 +0 - 13 
± - ZD -o.ii 

i 1^+0.28 
1 - iO -0.22 
1 14+0.23 
± - ± ^-0.21 


0.33 
0.43 
0.57 
0.52 


0.99 
0.99 
0.82 
0.88 



a mean interior slope of d log E/d log ~ —0.5). This in- 
ner profile matches well with the WL data in the region of 
overlap, at around ~ 200 kpc. The two SL modelling meth- 
ods cannot be distinguished in terms of the radial profile or 
the overall large-scale mass distribution, but differ mainly 
on fine scales principally owing to the inclusion of cluster 
members in the parametric method. 

Comparisons between other, different mass-modelling 
methods have been made before (e.g., Valls-Gabaud et al. 
2006, Coe et al. 2008, 2010, Donnarumma et al. 2010, 
Meneghetti et al. 2010). Here we find that the galaxy con- 
tribution is important to include in order to obtain ac- 
curate predictive power, as is the case in the parametric 
model (Zitrin et al. 2009b, Zitrin et al. 2010) which identi- 
fies and reproduces many multiple images, even if initially 
constrained only by a few obvious systems. On the other 
hand, being correlated to the initial light distribution, the 
parametric model may at times be less flexible and therefore 
over- sensitive to local luminous clumps or substructure, than 
the non-parametric model which does not make any prior as- 
sumptions about the input distribution of mass (Liesenborgs 
et al. 2006). 

We have made use of the reproduction ability of our 
parametric-model and the remarkable multiply-lensed ring- 
like system next to the cD galaxy, to uniquely constrain 
its projected mass and profile in the inner region ~ 1 — 5" 
(~ 4 — 22 kpc). We have also found that the low M/ Lb of 
~ 6 (M/L)q in this region can be fully accounted for by 
stars, similar to cD galaxies in other well-known clusters. 
The effect of baryons on halo density profiles is still unclear 
and may be related to "overcooling" claimed in studies of 
this effect (Barkana & Loeb 2010, Duffy et al. 2010, Mead 
et al. 2010). Further charactisation of cD galaxies from de- 
tailed lensing work will shed more light on this still poorly 
understood class of objects. 

We have found that Abell 1703 lies above the stan- 
dard c-M relation (Figure p"7|, similar to several other well- 
known clusters for which detailed lensing-based mass profiles 
have been constructed, adding to the claimed tension with 
the standard ACDM model (e.g., Broadhurst et al. 2008, 
Umetsu et al. 2010; see also Sadeh & Rephaeli 2008). Still, 
the overall level of systematic uncertainties may be too large 
to allow a definite conclusion regarding a clear inconsistency 
with ACDM predictions. To further explore this apparent 
discrepancy a substantial multi-cycle Hubble program has 
been established, based on an X-ray selected sample of re- 
laxed clusters so that no lensing bias is present in their se- 
lection (the CLASH program; see also §1), for which the 
SL techniques applied here will be of great value for deriv- 




Figure 17. Chi-squared confidence levels (grayscale; 68.3%, 
95.4%, and 99.7% confidence levels) of the joint SL±WL data fits 
to an NFW profile of Abell 1703, presented on the c-M plane. 
Overplotted are the expected c-M relations and their la uncer- 
tainties, presented in Duffy et al. (2008) with WMAP5 parame- 
ters, scaled to z c = 0.28. The cyan shaded curve corresponds to 
the full Duffy et al. (2008) sample, and the blue shaded curve 
corresponds to their relaxed-halo sample. As can be seen Abell 
1703 lies above the standard c-M relation. 



ing a statistically large and unbiased measurement of the 
equilibrium mass profiles of galaxy clusters. 
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